Connecting the vulcanization transition to percolation 

Weiqun Peng 1 , Paul M. Goldbart 1 and Alan J. McKane 1,2 
1 Department of Physics, University of Illinois at Urbana- Champaign, 1110 West Green Street, Urbana, Illinois 61801, U.S.A. 
2 Department of Theoretical Physics, University of Manchester, Manchester M13 9PL, UK 

The vulcanization transition is addressed via a minimal replica-field-theoretic model. The appro- 
priate long-wave-length behavior of the two- and three-point vertex functions is considered diagram- 
matically, to all orders in perturbation theory, and identified with the corresponding quantities in the 
Houghton-Reeve- Wallace field-theoretic approach to the percolation critical phenomenon. Hence, it 
is shown that percolation theory correctly captures the critical phenomenology of the vulcanization 
transition associated with the liquid and critical states. 
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I. INTRODUCTION AND OVERVIEW 

The vulcanization transition (henceforth denoted VT) is the generic equilibrium phase transition from a parent 
liquid state of matter to an amorphous solid state, driven by the imposition of a sufficient density of permanent 
random constraints between the constituents of the liquid. In the most common setting of the VT, the constituents 
of the liquid are macromolecules, the locations of which provide the annealed random (i.e. thermally equilibrating) 
variables. The constraints are commonly provided by covalent chemical bonds (i.e. cross-links), and impart quenched 
randomness upon the system. Over the past few years, a rather detailed description of the VT has been developed, 
ranging from a mean-field theory of the emergent amorphous solid state (including its structural and elastic 
properties ||, and its stability for reviews, see Refs. ||||) to the critical properties of the VT itself [Q. 

The present Paper aims to extend the description of the critical properties of the VT by exploring its relationship 
with the percolation transition (henceforth denoted PT; for a review, see Ref. Q). This relationship, which has long 
been anticipated on physical grounds ]^Jlfj| | , has recently found support both at the mean- field level , and beyond, 
via a renormalization-group approach ^^Specifically, it was recently shown that the order-parameter correlator near 
the VT (which probes for relative localization of particles) and its physical analog in PT (viz., the connectedness 
function) are governed by the same critical exponents, at least to first order in an expansion about the upper critical 
dimension six Q. 

The central result of the present Paper is the explicit reduction of certain basic critical properties of the VT to 
equivalent basic critical properties of the PT pT| . As we shall see, this reduction can be accomplished via an exact 
diagrammatic analysis of the complete perturbative expansion of the appropriate vertex functions of the VT. These 
are shown to furnish, in the replica limit, precisely the field-theoretic formulation of the PT due to Houghton, Reeve 
and Wallace (which henceforth we shall refer to as HRW). Hence, we establish that the critical properties of the 
VT and the PT are identical, not just to first order but to all orders in the departure of the spatial dimension d from 
the upper critical dimension. 

It is worth observing that the VT and the PT do, nevertheless, represent distinct physical phenomena. This is 
exemplified, e.g., by the amorphous solid state that emerges at the VT, which does not have an evident counterpart 
in PT. Another point of distinction is revealed by the role of fluctuations in low-dimensional systems, which are 
expected to have qualitatively different impacts on the states emerging at the VT and the PT jl3|. Yet another 
point of distinction concerns the nature of the degrees of freedom involved in the description of the VT and the PT. 
The former arises in systems having both quenched and equilibrating randomness, whereas the latter takes place in 
systems involving just one type of randomness (typically taken to be the quenched randomness); see Ref. [Jl3|. 

After completion of the present work we learned of the elegant work of Janssen and Stenull jl^], conducted inde- 
pendently of and simultaneously with the present work, which builds on earlier work on random resistor networks and 
percolation to arrive at, inter alia, essentially the same results as those contained in the present paper via a related 
approach. 

II. MINIMAL MODEL OF THE VULCANIZATION TRANSITION 

Our analysis is based upon a minimal model of the VT that accounts for thermal fluctuations in the positions 
of the constituents of the parent liquid, short-range repulsions between these constituents, and permanent random 
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constraints (e.g. resulting from cross-linking) which explicitly reduce the collection of configurations accessible to 
the constituents. This minimal model yields a rich mean-field picture of the structure and elastic response of the 
amorphous solid state, the former aspect having been verified by the computer simulations of Barsky and Plischke [ fL5| . 

The minimal model for the VT can be built (in the spirit of the Landau- Wilson scheme for continuous phase transi- 
tions) on the general basis of symmetry considerations and a gradient expansion 0]. The appropriate order parameter 
fl, whose expectation value detects the emergence of the amorphous solid state, has been discussed elsewhere [ |l6| . 
The quenched random constraints are accounted for via the replica technique, which incorporates the Deam-Edwards 
model fin\ for their statistics (viz. the statistics of the random constraints are determined by the instantaneous cor- 
relations of the unconstrained system) . The additional replica associated with the Deam-Edwards model leads to a 
situation in which one considers the n — > limit of a system of not of n but of n + 1 replicas. 

The resulting Landau- Wilson effective Hamiltonian takes the form of a cubic field theory involving the field fl, the 
argument of this field lying in (n + l)-fold replicated d-dimensional space: 

m = y^\ E (r + k ■ k^ Q (k)\ 2 + £ 0(*i) Sl(h) • 

feeHRS fei,fe 2 ,fe3€HRS 

The free-energy density / is (up to uninteresting factors that we shall ignore) related to this Hamiltonian via / oc 
lim„^o n _lm {/ DO. exp(— S)}; the functional integral is taken over the independent components of fl that feature in 
S. In S, the quantity ro is the VT control parameter which, near the VT, is linearly related to the density of random 
constraints. (To ease comparison with the HRW field theory of the PT, the coefficients and fields featuring in S are 
not defined exactly as they have been in our earlier works JTsj].) The symbol k denotes the replicated wave vector 
{k°, k 1 , . . . , k n }; the extended scalar product k ■ c is denoted k° • c° + k 1 • c 1 + • • • + k" • c™. The specification HRS 
arises from the following considerations. Consider the space of replicated wave vectors k. We decompose this space 
into three disjoint sets: (i) the higher replica sector HRS, which consists of those k containing at least two nonzero 
component-vectors k Q ; (ii) the one replica sector IRS, which consists of those k containing exactly one nonzero 
component-vector k"; and (iii) the zero replica sector ORS, which consists of the vector k — 0. This decomposition 
is illustrated schematically in Fig. [I] for the case of two replicas. It is especially straightforward to visualize this 
decomposition if the volume of the system is kept finite (and periodic boundary conditions are imposed) so that 
replicated plane waves with discrete, equally spaced, replicated wave vectors provide the natural complete set of 
functions. 




FIG. 1. Decomposition of the space of replicated wave vectors. Off-axis wave vectors lie in the HRS; on-axis (but off-origin) 
wave vectors lie in the IRS; the wave vector at the origin is the ORS. 

The symbol X^I-gHRS indicates a summation over the replicated wave vector k, subject to the restriction that k lies 
in HRS. This condition on the summation over k is essential. Physically, it reflects the fact that no macroscopically 
inhomogeneous modes (such as crystalline modes) order or fluctuate critically in the vicinity of the VT, such modes 
being stabilized by the excluded-volume interactions. Mathematically, this condition reduces the symmetry of the 
theory from one that contains the rotation group in replicated (n + l)<i-dimensional space to one that contains only 
rotations within each individual replica, along with the permutation of the replicas. 



2 



III. DEMONSTRATING THE EQUIVALENCE OF THE CRITICAL PROPERTIES OF THE 
VULCANIZATION AND PERCOLATION TRANSITIONS 



A. Overall strategy 

We now explain the strategy that we shall use to relate the VT and the PT. We shall focus on the replica limit 

(2) " f3l " 

of the long-wave-length behavior of the two- and three-point vertex functions, T n (k) and r„ (fci,^) in the VT 
field theory. The physical significance T n 2 \k) as a probe of connectedness has been elucidated in Ref. @. Now, the 
symmetry of the VT field theory dictates that the only suitably invariant term quadratic in the wave vector k is k ■ k. 

(2) " 

Thus, in a long- wave- length expansion for T„ (k), we have 

™ = r ?'<«> + whr d i (sf ■ sO tfwL k ■ k + ■ • • (3 - la > 

= A n + B n k ■ k H . (3.1b) 

As the upper critical dimension for the VT is six, and this is the dimension about which one may imagine expanding, 

(2) (3°) 

general renormalizability considerations demand that just these two vertex functions (r„ and P„ ) contain the 

primitive divergences, and do so via the constants A n , B n and C„ = F„ (6, 6) (see Ref. [p^|). 

Having identified the quantities central to a renormalization-group analysis of the VT, we shall establish that these 
quantities are identical, in the replica limit, to the corresponding quantities in percolation theory. To do this we shall 
make use of a convenient representation of the critical properties of percolation theory, viz., the HRW field theory 
representation (H. So that we know what we need to make contact with, we pause to give a brief account of this 
HRW representation, the Landau- Wilson Hamiltonian for which is given by 

H = |^|i(V0) 2 -i(V^) 2 + iro(0 2 -^) + |(0 + ^) 3 }, (3.2) 

where <fi is an ordinary field but ip is a ghost field. As HRW have shown, provided one enforces the rule that only 
graphs that are connected by <p-lines are included, the two- and three-point <f> vertex functions are identical (order by 
order in perturbation theory in the coupling constant g) to those of the one-state (i.e. percolation) limit of the Potts 
model. We mention, in passing, that this HRW representation consists of fields residing on c?-dimensional space, and 
does not necessitate the taking of a replica (or Potts) limit. However, it does require the additional rule by which 
certain diagrams are excluded by hand. 

Our strategy is as follows. Consider the standard Feynman diagram expansion for the two- and three-point vertex 



functions of the VT field theory in powers of the coupling constant g in the Hamiltonian (2.1). 

(i) To deal with the constraint that the internal wave vectors in the resulting diagrams reside in the HRS, we relax 
this constraint on summations over internal wave vectors but compensate for this by making appropriate subtractions 
of terms. 

(ii) Next, we observe that all diagrams for the two- and three-point vertex functions can be organized into two 
categories: those in which there is at least one route between every pair of external points via propagators having 
unconstrained wave vectors (which we call freely- connected diagrams); and the remaining diagrams, in which there is 
at least one pair of external points between which no paths exist consisting solely of propagators having unconstrained 
wave vectors (which we call freely-unconnected diagrams). Having made this categorization, we show that the appro- 
priate version of wave-vector conservation renders the freely-unconnected diagrams negligible in the thermodynamic 
limit, leaving us with a representation that is already reminiscent on the HRW approach. 

(iii) At this stage we have reduced the construction of the two- and three-point vertex functions to the computation 
of freely-connected diagrams only. Next, via a straightforward combinatorial analysis, we show that, in the replica 
limit, only a small class of diagrams survive. 

(iv) Finally, we explain how, again in the replica limit, the values of the remaining diagrams are precisely those 
occurring in the HRW prescription for percolation. 

We now set about implementing this strategy. 



B. Relaxing the constraint to higher replica sector wave vectors 



In the VT field theory the internal wave vectors occurring in the Feynman diagrams are constrained to lie in the 
HRS. In order to perform summations over these wave vectors, it is convenient to work with the continuum of wave 
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vectors (i.e. to take the thermodynamic limit) rather than the discrete lattice of them. In order to be able to take 
this limit, we re-express summations over HRS wave vectors in terms of unconstrained summations over (n + l)d- 
dimensional wave vectors, together with further unconstrained summations over (i-dimensional wave vectors (and 
certain trivial additional terms). To do this, we note that for a generic function F(k) we have 



E F(k) = £ i - ( £ 5>-, q n fc».«0 - 5 ko I m 

fceHRS k \ a=0q#O /3(jta) 



= £ I 1 - n 

Q=0 q /3(#a) 



^F(i)-^^F(0 ) ..„0,q« ) ) . 

£ a=0 q" 




(3.3a) 

(3.3b) 
(3.3c) 



which effects the re-expression of the summations just described. Note that, as it always comes with the factor n, 
the Sjj. q term will vanish in the replica limit, and can therefore be safely ignored. We shall refer to the wave vectors 
included in the term ^ Q as lower replica sector (LRS) wave vectors. Via these steps one can relax a constrained 
summation over HRS wave vectors, instead freely summing over all replicated wave vectors, provided one compensates 
by augmenting the summand with the factor 



i-EE^ n <w 

a=0 q p&a) 



(3.4) 



How does this constraint relaxation manifest itself in the setting of Feynman diagram computations? One simply 
augments every internal propagator V n+1 Go(k) with a factor (EOI): 



Go(fe) = 



1 



ro + k ■ k 



ro + k ■ k 
1 



r a + k-k r + k ■ k a=Q q 



n <w- 



(3.5a) 
(3.5b) 



How this decomposition is expressed diagrammatically is shown in Fig. ^. 



H = 







FIG. 2. Decomposition of an HRS propagator (indicated by H) into an unconstrained propagator (the unadorned line), less 
an LRS propagator (indicated by L). 



In this manner, each HRS internal line in a Feynman diagram can be decomposed into an unconstrained internal 
line, less an LRS internal line. Thus, the various vertex functions can be expressed in terms of Feynman diagrams 
composed of unconstrained lines together with as LRS internal lines. Note, for future reference, that physically 
meaningful vertex functions have external wave vectors in the HRS. 

We illustrate this decomposition for the case of a simple diagram in Fig. ||. More generally, we arrive at the following 
modified Feynman rules for the VT field theory: 

(i) Write down all diagrams arising from the original theory. In these diagrams all wave vectors are constrained to 
the HRS. 
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(ii) Replace each diagram with the collection of diagrams obtained by allowing each internal line to carry cither 
an unconstrained replicated wave vector or a LRS wave vector. (Thus, a diagram with L internal lines spawns a 
total of 2 L diagrams.) Identically- valued diagrams can be represented by a single diagram, together with a suitable 
combinatorial factor (see, e.g., the factor of 2 in Fig |J). 

(iii) Provide a factor of —1 for each LRS internal line. 

At this stage we observe that the combinatorics of our diagrammatic expansion coincides with those of the HRW 
expansion, provided one identifies the internal unconstrained and LRS lines of the VT theory with, respectively, the 
corresponding internal <p- and ?/>-lines of the HRW representation. 




FIG. 3. Decomposition of the one-loop diagram for the two-point vertex function. On the right-hand side of this equation, 
the first and second diagrams are freely-connected diagrams, and the third is a freely-unconnected diagram. 

We have, however, yet to show that the diagrams removed by hand in HRW can be safely omitted from the VT 
theory, and that the numerical values of the (replica limits of) the VT diagrams are identical to those of the HRW 
diagrams. We shall establish these facts in the following subsections. 



C. Elimination of freely- unconnected diagrams 

We remind the reader that in the HRW theory for the two- and three-point <^-field (i.e. the physical) vertex functions, 
one is instructed to remove, by hand, those diagrams in which there is at least one pair of external points between 
which no paths exist consisting solely of 0-field propagators. The corresponding diagrams in the VT theory are those 
in which there is at least one pair of external points between which no paths exist consisting solely of propagators 
having unconstrained wave vectors, i.e., freely- unconnected diagrams. We now show that these freely- unconnected 
diagrams of the VT theory automatically vanish in the thermodynamic limit. 
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FIG. 4. Schematic illustration of a piece of a diagram in the VT field theory obtained by cutting LRS lines in a 
freely-unconnected diagram. Note that the wave vector k flowing through the external point lies in the HRS. The shaded 
circle represents any way to connect the exhibited lines using the cubic interaction vertex, unconstrained propagators and LRS 
propagators. 

To do this, consider a generic VT theory diagram for the two- or three-point vertex functions. Observe that freely- 
unconnected diagrams have the following property: as there exists a pair of external points not connected by a path 
of unconstrained internal lines, there must exist at least one scheme of cutting solely LRS internal lines that causes 
the diagram to separate into disconnected pieces with the external points shared amongst the pieces. As we are 
considering only two- and three-point vertex functions, at least one of these pieces involves only a single HRS external 
point, along with a number of cut LRS lines. A schematic illustration of such a piece is shown in Fig. 

Let us examine the consequences of applying wave-vector conservation to this piece, noting that the wave vector 
k flowing in through the external point must lie in the HRS, whereas the wave vectors flowing out through the 
remaining (i.e. cut) lines lie in the LRS. Now, according to the VT field theory, wave-vector conservation requires 
that the incoming HRS wave vector k be equal, replica by replica, to the sum of the outgoing (m = 2,3, . . .) LRS 
wave vectors flowing in a given replica, i.e., that 

m 

k a =J2 sa ' a3 Pj (for a = 0,1,2,..., n) (3.6) 

3=1 

where cvi,a;2, . . . indicate the replicas through which wave vectors Pi,P2, ■ ■ ■ flow. As a consequence, because the 
incoming wave vector lies in the HRS, the outgoing LRS wave vectors must flow out through more than one replica. 
This is the key observation, as the following example, depicted in Fig. [I| reveals. Here, there are six outgoing LRS 
lines, two with wave vectors flowing in replica 0, one with wave vector flowing in replica 1, and three with wave 
vectors flowing in replica 2. For replica 0, wave-vector conservation reads k° = pi + p 2 , so that, e.g., pi determines 
P2- Similarly, for replica 1, wave- vector conservation reads k 1 = P3, so that P3 is determined. More generally, 
as this special case exemplifies, the number of independent outgoing LRS wave vectors is reduced by at least two 
(rather than one that total wave- vector conservation demands) simply because of the fact that the outgoing LRS wave 
vectors must flow out through more than one replica. This, in turn, means that in the uncut diagram there are fewer 
independent wave vectors to be summed over than the number of loop wave vectors suggested by simple topological 
counting. As a result, additional denominators of V n+1 remain, even after the summation over independent wave 
vectors in the uncut diagram are replaced by their thermodynamic-limit integrals, which renders the corresponding 
freely-unconnected diagram negligible. 

As a concrete example of the foregoing argument, we compute the third diagram on the right-hand side of the 
equation depicted in Fig. [| In this diagram, both of the internal lines lie in the LRS, the diagram does not (by simple 
wave-vector conservation) contribute unless the external wave vector k has nonzero d-vector components in precisely 
in two replicas (e.g. replicas one and two). In this case, the diagram makes the contribution 

2V n+1 G (k 1 ) V" +1 G (k 2 ) ( 5 r 2(n+1 ») 2 = 2 5 2 V- 2 (' l+1 )Go(k 1 )G (k 2 ). (3.7) 

On the right-hand side, one denominator of V n+l will combine with the Kronecker 5-function to maintain overall 
wave- vector conservation (via a Dirac 5-function in the thermodynamic limit); the other V n+1 denominator (which, in 
usual cases, would combine with the summation over a loop wave vector to produce an integral) makes this diagram 
vanish. This special case exemplifies the general emergence, in the VT setting, of the central aspect of the HRW 
formulation, viz., the removal of the (^-unconnected diagrams. 

D. Replica sums and their decomposition in the replica limit 

Now that we have demonstrated that only the freely-connected VT field theory diagrams contribute, we make a 
closer examination of these diagrams for the relevant cases of the two- and three-point vertex functions. We begin 
by noting that each diagram exhibits the full symmetry of the VT field theory, viz., invariance under separate d- 
dimensional rotations in each replica and permutations of the replicas. Therefore, the small-wave- vector expansion of 
Tn^ (fc) given in Eq. ( |3.1b| ) remains valid, diagram by diagram. 

Now, the computation of any contributing diagram involves summations over a number of independent LRS (but 
otherwise unconstrained) wave vectors as well as summations over the replicas through which these wave vectors flow. 
These latter summations over replicas can be decomposed as follows: 



G 



^ t Fa 1 .a 2 ,---,ai — ^ ] + ^ ] + ^ ] +•■■ ^ ] Fa 1 ,a 2 ,...,ai, (3-8) 

ai,02,...,Ql \ °1' Q 2 °>l c»l,Q2. --.°l a I ,a 3 ,...,a| a ll Q 2 >- ■.»( / 

\ all equal two distinct three distinct all distinct / 

where F ai a2> ... <ai is a generic function of the I replica indices. Said equivalently, the summation can be decomposed 
into: terms in which all wave vectors flow through a common replica; those in which the wave vectors flow through 
two distinct replicas; etc. 



Now let us make use of this decomposition. For A n and C n in Eq. ( 3.1t ) the external wave vector is zero [p.9| , 
and therefore the summand F ai ,a a ,...,ai ^ s invariant under permutations of the replicas J2(i| ]. Thus, in the first term of 
the decomposition, F is constant [i.e. independent of the (common) value of the replica arguments], and hence this 
term contributes (n + l)i*b,o....,o- In the replica limit, this becomes Fb,o,. ..,o- As for the second term, let us further 
decompose it into partitionings of the set of replica indices into two subsets, the replica indices in each subset having 
a common value. In each such partitioning F is constant, and thus each partitioning contributes (n + l)nF, which 
vanishes in the replica limit. By continuing with this decomposit ion tactic via tri-partitioning, tetra-partitioning, etc., 



we establish that all of the terms on the right-hand side of Eq. (3.8) except the first vanish in the replica limit. 

We next consider the coefficient B n . As mentioned above, symmetry considerations dictate that each diagram 
contributing to the two-point vertex function has the small-wave-vector expansion 

+ B ( f^k-k + --- , (3.9) 

where A^ 111 ^ and B^'^ are the contributions to A n and B n from the diagram in question. We exploit the (larger 
than mandated) (n + l)<i-dimensional rotational invariance of the terms retained in this small-wave-vector expansion 
by choosing k to be rotated into a single replica: {k, 0, . . . , 0}. (Although it has, until this stage, been vital to ensure 
that k lies in the HRS, e.g., in order to eliminate the freely- unconnected diagrams, one is now at liberty to ignore 
this requirement.) Repeating the tactic just used for the analysis of A n and C„, with the slight elaboration needed to 
accommodate the fact that the (suitably rotated) external wave vector {k, 0, . . . , 0} breaks the permutation symmetry 
group from V n +i down to V n . In this way, we see that the only contributions that survive the replica limit are from 
the all-equal partition and, furthermore, from the case in which all of the independent LRS wave vectors lie in replica 



E. Feynman integrals and their reduction to HRW integrals in the replica limit 



Having shown that the topology and combinatorics of the VT field theory diagrams for A n , B n and C„ coincide with 
those of the HRW field theory diagrams, the task that remains is to show that, for every diagram contributing to these 
coefficients, the actual value of the corresponding Feynman integral reduces, in the replica limit, to the appropriate 
HRW value. That this is so can most straightforwardly be seen by employing the Schwinger representation [ElJ of the 
powers of the propagator, viz., 



(*-!)! 
(ro + k ■ ky 



a s-l e -a(r + k-k) = / a s-l c -ar„ J| g-^-k" ^ g = ^ ^ ^ (3 ^ 



a=0 



where the Schwinger parameter a ranges between and oo. Observe that Eq. ( 3.10| ) presents the propagator in a 
form that is very conveniently factorized on the replica index. 




FIG. 5. A two-loop diagram used to exemplify the decoupling of the replicas using the Schwinger representation. 
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Let us begin with a concrete example. In the thermodynamic limit, the diagram depicted in Fig. |B] contributes to 
A n a term proportional to 



g 4 J S n+1 ^ d k d d p G Q (k) 2 G (p) 2 G (fc + p) 2 

= g 4 [ d d k° d d p [ CTl(T2CT3 e-( CT1 + CT2 +^)'-°e-( CT1 + CT3 ) k0 - k0 e-^ +CT3 ) p - p e- 2CT3k0 - p f[ f d d k a c -^i+^)k°-k° 

J J o~i&2&3 a=l 

= g 4 fd d k°d d p f CTl a 2 a3e-( CTl + ff2 + ff3 ) ro e-( CTl + ff3 ) k ° k0 c-^+ ff3 )P Pc- 2ff3k0 P ( f d d k G -^+ a ^X 

g 4 [ d d kd d p { o- 1 o-2(T3e- ( - CTl+CT2+ff3)ro e-( ffl+CT3 ) k - k e-^ 2+ff3 ) p - p e- 2ff3k - p 
= 5 4 y d d fc d d p G (k) 2 G (p) 2 G (k + p) 2 . (3.11) 

This limiting value is precisely that occurring for the corresponding diagram in the HRW field theory for the PT. 

The tactic that we have just employed, viz., the use of the Schwinger representation to decouple the replicas from 
one another, provides easy and explicit access to the replica limit and, hence, to the precise correspondence with the 
HRW prescription. It can straightforwardly be invoked not only for all diagrams that contribute to the coefficients 
A n and (by the same procedure) G„, but also for the coefficient B n . 

When considering diagrams contributing to A n and G n , we saw that what survived were terms in which all internal 
LRS wave vectors flowed in a common (but otherwise arbitrary) replica. Now, as we consider _B„, there is a slight 
complication arising from the presence of an external wave vector, which spoils the full Vn+i permutation symmetry. 
However, this external wave vector has been chosen to lie in replica zero and, as we have shown above, the only 
surviving contribution is the one in which all internal LRS wave vectors also flow in replica zero. Then, via the 
Schwinger representation of the propagators, and via factorization on the replica indices, we see that the Feynman 
integral, in the replica limit, is identical to that in the HRW approach. Hence, the VT presents not only the the 
same coefficients Aq and Go of the two- and three-point vertex functions as does the HRW representation, but also 
the same coefficient Bq. 



IV. CONCLUDING REMARKS 



Let us summarize what is presented in this Paper. We have addressed the vulcanization transition via a minimal 
field-theoretic model. This model is built from an order parameter whose argument is the (n + l)-fold replication of 
ordinary d-dimensional space. [The structure of this theory should be contrasted with that of more familiar replica field 
theories, in which it is the (internal) components of the field that are replicated rather than the (external) argument] 
We have considered appropriate long- wave-length aspects of the two- and three-point vertex functions for this model, 
to all orders in perturbation theory in the cubic nonlinearity. Via a detailed analysis of the diagrammatic expansion 
for these quantities, we have found that, in the replica limit, these vulcanization-theory vertex functions precisely 
coincide with the corresponding vertex functions of a certain field-theoretic representation (due to Houghton, Reeve 
and Wallace) of the percolation transition. Hence, percolation theory correctly captures the critical phenomenology 
the liquid and critical states of vulcanized matter. 
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